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STEADY-STATE  MODEL  OF  A FLAT  LASER-DRIVEN  TARGET 


In  an  earlier  paper^  exact  transmission  coefficients  were 
calculated  for  intense  light  incident  on  a plasma  slab  in  which  ions 
were  frozen.  In  this  paper  a steady-state  model  of  a plasma  slab, 
accelerated  by  its  interaction  with  laser  radiation,  is  treated  by 
1-D  steady-flow  hydrodynamic  equations  in  an  accelerated  frame  of 
reference.  The  cold  unablated  fluid,  the  ablation  layer,  both 
classical  and  flux-limited  hot  conduction  regions,  the  critical 
surface,  and  the  underdense  blowoff  are  all  considered.  A global 
description  determines  the  temperature,  density,  velocity,  and 
boundaries  as  well.  Approximate  analytic  solutions  are  given  In  each 
of  the  plasma  regions. 

The  ablation  layer,  containing  a steep  density  gradient  separating 
cold  dense  fluid  from  hot  low  density  plasma,  moves  at  the  front  of 
the  thermal  wave  and  is  accelerated  by  the  rocket  reaction  to  the 
ablation.  Nevertheless  J.  Boris  has  demonstrated  through  time- 
dependent  numerical  simulations  that  the  temperature  profile  at  the 
ablation  layer  is  steady  in  an  accelerated  reference  frame,  and  has 
provided  a simplified  analytic  model  of  the  steady  temperature  profile 
at  the  ablation  layer. ^ 


The  quasisteady-state  model  presented  here  accounts  self 
consistently  for  the  (slow)  increase  In  acceleration  of  the  slab 
owing  to  the  diminishing  mass  of  cold  fluid  being  accelerated. 

While  a slab  geometry  reveals  most  of  the  important  physics,  and 
can  be  used  to  consider  flat  target  experiments  and  instability 
at  the  ablation  layer,  the  model  readily  transforms  to  a 
spherical  geometry.  Max  et  al?  have  found  analytic  steady-state 


Note:  Manuscript  submitted  July  18,  1977. 


solutions  for  plasma  flow  In  the  hot  conduction  region  of  spherical 
laser  targets,  and  Gitomer  et  al.**  have  treated  spherical  flow  with  an 
artificial  mass  source. 

The  model  is  formulated  in  terms  of  time-independent  solutions 
of  the  continuity,  momentum,  and  energy  equations. 


0 


_ ^ 

" at 


^ A 1+21 

0 = ^ (DV)  = - ^ (pv2  + P)  + pg  . (-2— I)5(x  . , 

0 = II  = - ^ (Pv  + ev  + q)  + pgv  + I^8(x  - x^). 


(1) 


Here  v is  the  velocity  of  a volume  element  of  mass  density  p, 
temperature  (in  energy  units)  T,  pressure  P = pT/m^,  ion  mass  m^,  and 
local  energy  density  6 = (5P  + pv^)/2;  absorbed  (reflected)  laser  flux 
is  1,  (V. 

The  laser  energy  and  momentum  is  regarded  as  deposited  at  the 
critical  surface  x = x^.  The  momentum  deposition  term  causes  gradient 
steepening  and  density  shelves  at  the  critical  surface  that  have  been 
observed  experimentally,^  and  calculated  theoretically®  by  other  means. 


The  effective  gravity  g must  be  included  in  the  energy  equation 
as  the  product  of  force  density  pg  with  velocity  v.  It  cannot  be 
omitted,  as  in  Ref.  7. 

The  heat  flux  q at  any  point  in  the  hot  conduction  region  is 
considered  to  be  the  lesser  in  magnitude  of  the  classical  flux 
%lass  " and  limited  flux  = - (,p(T/mj^)^'^^T'/l  t'|  . A 

prime  denotes  d/dx,  KT^^^  is  conductivity,  and  the  flux-limit  parameter 
lies  in  the  author-dependent  range^  0.5  ^ I ^ 60.  An  upper  bound, 

•t  < 3,  for  this  model  is  derived  below. 

The  continuity  equation  implies  constant  momentum  density  pv  = p^v^. 
Fluid  quantities  evaluated  at  the  origin,  chosen  at  the  surface  of 


2 


maximum  density,  are  denoted  by  a subscript  zero.  Then  the  momentum 
and  integrated  energy  equations  in  the  overdense  fluid  simplify  to 
two  first  order  equations  in  p and  T 

p'  = p(T'  - m^g)/rm^v^(p^/p)^  - T]  , 

q(T,T',p)  = - ^(PjjV^/m^)r5(T  - T^)  + mj^v|(p2/p2  - 1)  - 2m^gx] 

that  can  be  solved  by  a Runge-Kutta  code  in  classical-conduction  and 
flux-limited  regions.  Solutions  near  the  ablation  layer  are  shown 
in  Fig.  1. 

c /2 

The  flux  ~ T^  T'  in  the  cold  fluid  rises  monotonically  towards 
the  heat  source,  but  the  temperature  gradient  T',  while  positive 
definite,  is  not  necessarily  monotonic.  At  maximum  density,  T^  = m^g; 
conduction  is  enhanced  by  increased  acceler&tion. 

If  the  dimensionless  gravity  ^ m and  Mach  number 

M s (m^v^/T)^  are  both  « 1 in  the  cold  fluid,  then  approximate 
analytic  solutions  of  (2)  are 


T « T^(l  - 6r/25  + 2m^gx/5T^)  , 

p » p^(l  + 3m^gx/5T^)  , (-T^/m^^g  < x < 0)  . 


O) 


Indeed  the  laminar  ablation  rate  is  generally  much  less  than  the  local 
sound  speed  at  the  ablation  layer  in  ablative  implosions,  though  not 
in  explod ing-pusher  targets. 

A sonic  point  (M  = 1)  can  occur  in  the  cold  fluid  at 


T 

0 

2 "\p(x3); 

3 


If  M « 1,  and  T « 1,  Chen  an  approximate  analytic  solution  of 

(2)  in  the  hot  classically  conducting  region,  valid  for  0 < x < T /m.g 

o 1- 

is 


4m 


.kt5« 


x(T)  W 


25  P V 
o o 


- 1)  +1  - 1)  + 5 (T^  - 1)], 


in  which  t 9 T/T  . The  analytic  solutions  (3)  and  (^)  for  cold  and 
o 

hot  conduction  regions  are  compared  in  Fig.  1 with  the  numerical 
solution  of  (2).  At  high  temperatures  (^)  can  be  inverted 


T « (25d  V - 2T  /5  , (T  » T ) . 

'■  ^“o  o i o ’ o 


At  sufficiently  high  temperature  and  low  density,  the  heat  flux 

may  saturate  in  the  overdense  region.  The  solutions  of  (2)  for  both 

classical  flux  and  saturated  flux  are  matched  where  q , = q,  . . 

class  iim 

In  a saturated-flux  region  (2)  implies® 


/c  j.  m2 


The  gravity  gives  the  Mach  number  a spatial  dependence  and  boost; 
Che  flow  is  assisted,  and  onset  of  flux  limiting  delayed. 

At  the  critical  surface,  conservation  of  mass  requires  s .v.  = d v , 
so  that  the  jumps  in  density  and  velocity  are  related  by 


v_^Ap  = - p_Av  . (6) 

The  notation  concerning  the  jump  in  any  variable  u is  defined  by 
Au  s u^  - u_  s u(x^  + 0)  - u(x^  - 0).  All  energy  deposited  at  the 
critical  surface  is  assumed  to  be  conducted  into  Che  much  more  massive 
and  colder  overdense  region,  so  that  q = 0 in  the  underdense  region,  and 
an  adiabatic  equation  of  state  applies  there.  Then  with  AT  = C and 
T(x^)  * Aq  = “ q(p  » T^)  > 0.  The  jump  conditions  from  (1)  are 
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Av  = (cp  - l)v  , (7) 


Ao  = - m (I  + 21  )/(T  - m v2ce)c  , 

^ a ^ C i * 

where  ® a T 1 + 2(1^  - Aq)/p_v^]^.  The  case  is  considered  for  which 
the  heat  conducted  away  from  does  not  exceed  the  laser  energy 
deposited  there,  so  that  © > 1.  Together  (6)  and  (7)  imply 

+ 21^  = (p_c/mj^)(T^  - m^v2(p)(l  - l/^)  , (8) 

Since  I +21  >0  and  o > 1,  then  ra.v^  < T . Thus  steady  flow  into 

the  critical  surface  is  subsonic  In  the  rest  frame  of  the  critical 
surface.®  Then  (5)  implies  the  flux- limit  parameter  must  be  less 
than  3 for  a flux-limited  region  to  exist  in  steady  state.  This 
upper  bound  is  comparable  with  recent  experimental  values  and 
theoretical  predictions.®’"^ 

Mass,  momentum  and  energy  conservation  across  the  critical 
surface  insufficiently  constrain  the  solution.®  A reasonable  added 
constraint  derived  in  Ref.  6 for  a related  problem,  and  independent 
of  ponderomotive  force,  is  p_  = Pj.rMf  " -tnCM^)  - l‘!/2(l  - M_)^.  The 
density  shelves  and  subsonic  flow  into  the  critical  surface  are 
features  of  other  theoretical  treatments®  as  well,  and  are  seen  here 
to  be  caused  exclusively  by  laser  momentum  deposition. 

2 

Ordinarily  a ponderomotive  force  term  proportional  to  pdf eI  /dx 
is  appended  to  the  momentum  equation  to  model  the  effects  of  the 
electric  field  E in  the  underdense  region,  imposing  spatially  periodic 
fluctuations  on  density  and  velocity  profiles.  In  this  model  the 
momentum  equation  is  averaged  over  many  wavelengths  to  iron  out  the 
effects  of  the  ponderomotive  force,  (pd|E|^/dx)^  = C,  and  leave  only 
the  gross  profiles  of  the  steady  flow  of  an  ideal,  adiabatic  gas  in 
the  underdense  region. 

The  exact  solutions  of  (1)  for  the  underdense,  adiabatic  gas 

are 
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(9) 


- 5(T^/mpri  - (v^/v)^^5]  = + 2g(x  - x^) 

P = 0+%/^  , T = T^Cv+/v)^^^ 

Cn  physical  grounds  T must  decrease  away  from  the  heat  source 
at  x^.  Then  (9)  implies  the  flow  out  of  the  critical  surface  must 
be  supersonic,  > 5T^/3, 

The  gravity  determines  the  boundaries  self  consistently.  Let  the 
total  plasma  mass  per  area,  m,  be  specified.  Since  momentum  density 

is  uniform,  during  a time  6ti  a thin  layer  of  mass  6m  = o v 6t  = p v *t 

L L * * 

is  effectively  transferred  from  the  left  boundary  x^  to  the  right 
boundary  x^^  at  higher  velocity  v^^.  The  change  in  momentum, 

6p  = 6m(Vj^  - must  be  compensated  by  an  acceleration  of  the  whole 

plasma  to  the  left,  since  the  hydrodynamic  forces  are  internal, 
causing  an  effective  gravity 

g = P.v_(Vj^  - v^)/m  . (IG) 

This  condition  determines  the  approximate  slab  boundaries,  and 
completes  the  global  description  of  the  laser-plasma  interaction. 

Figure  2 shows  the  self-consistent  profiles  and  boundaries  found  by 
specifying  I^,  I^,  m,  g,  p^,  and  T^  only. 

Since  the  jump  conditions  are  not  satisfied  at  the  slab  boundaries, 
the  steady  state  assumption  breaks  down  near  the  boundaries.  Equation 
(10)  is  good,  however,  as  long  as  the  widths  of  the  unsteady  slab 
ends  are  much  less  than  the  overall  width  of  the  slab. 

The  gravity  increases, 
dg/dt  = p_v_(dv|/dx  - dv^/dx)/2m  , 

and  can  only  be  considered  steady  during  an  interval  At  if 
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(11) 


At  « 2(Vj^  - v^)/(dv|/dx  - dv^/dx)  . 

If  the  cold  fluid  is  no  wider  than  about  T /m.g,  and  the  underdense 

o 1 

region  is  wider  than  about  Tg/®^8.  then  (3),  (9),  and  (11)  imply  a 
characteristic  time  scale  for  growth  of  the  acceleration  of 
The  dynamical  development  of  a slab  can  be  followed  over  longer  time 
periods  by  respecifying  gravity  and  intensity  if  desired,  and 
repeating  the  integration  of  (2). 

I wish  to  thank  J.  P,  Boris  for  valuable  comments  and  suggestion 
of  this  topic,  and  S,  E.  Bodner  and  D.  L,  Book  for  advice  and 
encouragement.  ' 
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Fig.  1 — Temperature,  density,  and  pressure  profiles  at  ablation  layer 
calculated  by  model  for  T^/mj  = 10^^  cm^/sec^,  v,,  = 2.2  X 10®  cm/sec, 
Po  = 0.3  g/cm®,  and  acceleration  g = 3 X 10^®  cm/sec^.  Dashed  curve 
is  analytic  temperature  approximation  (Eqs.  (3-4)). 


Fig.  2 — Density,  temperature,  and  Mach  number  profiles  of  global 
slab  model.  Specified  were  absorbed  flux  Ij,  = 10.2  TW/cm^,  re- 
flected flux  If  = 10.2  TW/cm^,  total  mass  m = 0.733  mg/cm^,  accel- 
eration g = 3 X 10^®  cm/sec2,  critical  density  = 4 X 10“^  g/cm^, 
and  critical  temperature  T./mj_=  4.50  X 10^“*  cm^/sec^;  coefficient 
of  conductivity  K = 10“ ^ (cgs).  Self  consistency  of  the 
global  model  required  critical  surface  at  Xg,  left  and  right  boundaries 
at  Xl  and  xj^,  and  upper  and  lower  shelf  densities  p_  and  p+  shown. 
No  saturation  of  heat  flux  occurs  here  for  flux-limit  parameter  C > 


T/rtij  (cm^/sec^) 
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